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HYPERBOLIC  CONSERVATION  LAWS  IN  VISCOELASTICITY 


J.  A.  Nohel1,  R.  C.  Rogers2,  and  A.  E.  Tzavaras3 


1.  Introduction. 


The  equations  of  unsteady  motion  for  nonlinear  elastic  materials  are  quasilinear  sys¬ 
tems  of  hyperbolic  equations  in  which  characteristic  speeds  are  not  constant.  Thus  weak 
initial  waves  are  amplified  and  smooth  solutions  generally  develop  singularities  in  finite 
time.  Particularly  interesting  situations  arise  when  this  destabilizing  mechanism  coex¬ 
ists  and  competes  with  dissipation;  a  simple  example  is  provided  by  the  quasilinear  wave 
equation  with  linear  first-order  damping.  A  more  subtle  dissipation  occurs  in  certain  vis¬ 
coelastic  materials  such  as  polymers,  suspensions  and  emulsions  which  have  memory;  i.e. 
the  stress  at  each  material  point  depends  not  only  on  the  present  value  of  the  deformation 
gradient  (and/or  velocity  gradient),  but  on  the  entire  temporal  history  of  motion.  These 
materials  exhibit  behavior  intermediate  between  that  of  an  elastic  solid  and  a  viscous  fluid. 
Generally,  the  memory  fades  with  time:  disturbances  which  occurred  in  the  distant  past 
have  less  influence  on  the  present  stress  than  those  which  occurred  in  the  recent  past. 

J 

Our  purpose  is  to  discuss  qualitative  properties  of  the  equations  which  model  unsteady 
motions  of  such  materials.  In  order  to  present  some  of  the  central  ideas  while  avoiding 
serious  technical  difficulties,  we  shall  restrict  the  discussion  to  a  particular  model  problem 
for  the  motion  of  a  one-dimensional  viscoleastic  material  with  fading  memory.  Many  of  the 

I 

qualitative  features  in  this  relatively  simple  situation  are  present  in  more  general  integral 
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and  differential  models  in  one  and  in  several  space  dimensions;  these  are  discussed  in  a 
forthcoming  monograph  [37]. 

Motivated  by  the  study  of  the  growth  and  decay  of  acceleration  waves  by  Coleman 
and  Gurtin  ([5],  [6])  and  by  experiments  ([34j,  [43]),  it  is  natural  to  require  that  reason¬ 
able  mathematical  models  for  unsteady  motion  of  viscoelastic  materials  should  possess 
qualitative  properties  which  include: 

(a)  Under  physically  reasonable  assumptions,  the  memory  should  induce  a  weak  dissipative 
mechanism.  As  a  consequence,  for  sufficiently  smooth  and  small  data  and  history  the 
equations  of  motion  should  have  globally  defined  smooth  solutions  which  decay  as  t  — ►  oo. 

(b)  By  contrast,  if  the  smooth  data  and  history  are  chosen  sufficiently  large,  the  classical 
solution  solutions  should  break  down  in  finite  time  and  exhibit  a  shock  structure. 

(c)  For  rough,  bounded  data,  the  equations  of  motion  should  possess  globally  defined  weak 
solutions. 

This  paper  is  organized  as  follows.  In  Section  2,  we  state  the  model  problem  and  re¬ 
lated,  well-understood  problems,  and  we  briefly  review  known  results  regarding  properties 
(a)  and  (b)  for  classical  solutions.  We  then  state  a  new  result  [32]  on  the  existence  of  a 
weak  solution  (in  the  class  of  bounded  measurable  functions)  which  verifies  property  (c) 
in  an  important  special  case.  In  Section  3,  we  sketch  the  proof  of  the  new  result  using 
the  methods  of  vanishing  viscosity  and  compensated  compactness;  further  details  can  be 
found  in  [32]. 


2.  Model  Problem  and  Summary  of  Results. 


1 

& 


We  study  the  model  Cauchy  problem  for  the  system  of  integrodifferential  equations 


wt  =  vx,  vt  =  ax,  x  6  R,  t  >  0, 
w(x,0)  =  tvo(x),  v(x,0 )  =  t>o(x),  x  £  R, 


(V£) 


where  the  function  a(x,t)  is  determined  by  the  history  of  w(x,-)  through  the  constitutive 


assumption 


r(x,t)  =  <p{w{x,t))  +  f  k(t  -  r)V’(u;(x,r))dr, 

Jo 


( CA ) 


WWW 
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which  generalizes  Boltzmann’s  law  for  linear  viscoelasticity  [3].  The  given  functions  <p(u>), 
i/>(ti’)  and  k(t)  are  assumed  to  be  smooth  and,  in  addition, 


<p'{ tv)  >0,  ui£R, 

so  that  the  structure  of  (VE)  is  hyperbolic.  The  system  (VE)  is  a  model  for  one-dimensional 
motion  of  a  viscoelastic  material  with  fading  memory.  The  functions  v[x,t),  w(x,t)  and 
cr(x,t)  stand  for  the  velocity,  deformation  gradient  and  stress,  respectively,  while  the  con¬ 
stitutive  assumption  (CA)  states  that  the  stress  is  a  particular  functional  of  the  history  of 
the  deformation  gradient;  in  (CA)  the  history  of  w(x,  •)  is  assumed  to  be  zero  for  t  <  0 
and  the  body  force  is  taken  to  be  zero.  Under  appropriate  assumptions  on  the  kernel  k, 
the  memory  induces  a  weak  dissipative  mechanism  which  competes  with  the  hyperbolic 
character  of  (VE).  In  the  sequel,  we  will  limit  our  discussion  to  the  Cauchy  problem. 
We  remark  that  when  x  is  restricted  to  a  finite  interval,  boundary  conditions  which  are 
compatible  with  the  data  wo,v0  must  be  adjoined  to  (VE). 

It  is  useful  to  consider  several  special  cases.  If  k  =  0,  (VE)  reduces  to  the  system  of 
conservation  laws 

f  =  vx,  vt  =  <p{w)x 

\  u;(x,0)  =  te0(z),  v(x,0)  =  v0(x),  '  ’ 

which  describes  the  motion  of  one-dimensional  elastic  materials.  The  solutions  of  (E) 
generally  develop  singularities  in  finite  time  no  matter  how  smooth  and  small  the  data 
w o,vo  are  taken  ((I9],[24],[l8j).  At  the  other  extreme,  if  one  formally  sets  k  =  —S',  where 
6  is  the  Dirac  mass  at  the  origin,  (VE)  reduces  to  the  system 


wt  =  vx ,  vt  =  il>[w)xt  +  <p{w)x, 

which  is  parabolic  and  which  possesses  globally  defined  smooth  solutions  even  for  the  large 
data  u/o,  t>o,  whenever  x l>  is  smooth  and  t/’'(- )  >  0  ([l],[2lj).  The  theory  of  (VE)  lies  between 
these  two  extremes.  In  the  special  case  0  =  <j>,  (VE)  reduces  to  the  system 

vt  =  <p(w)x  +  fQ  k(t  -  r)p(w(x,r))xdr,  x  €  R,  t  >  0, 

(2.1) 

,  te(x,0)  =  w0(x),  v(x,0)  =  v0(x),  x  €  R, 
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which,  as  we  shall  see,  exhibits  a  close  resemblance  to  the  theory  of  the  system 


u>t  =  vx 

vt  =  <p{w)x  -  v. 


(FE) 


(FE)  describes  the  motion  of  a  one-dimensional  frictionally  damped  elastic  material  and 
enjoys  properties  (a),  (b)  and  (c)  which  are  established  in  [29,38,13]  respectively.  The 
similarity  between  (2.1)  and  (FE)  is  revealed  by  the  following  idea  of  MacCamy  [22].  Let 
r(t)  be  the  resolvent  kernel  associated  with  k;  i.e.,  r  is  the  solution  of  the  linear  Volterra 
equation 


r(t)  +  f  k(t  -  r)r(r)dr  =  k(t),  t  >  0. 

Jo 

Convolving  (2.1)2  with  r(t),  a  simple  calculation  yields 

[  k(t  -  t)<p(w(x,t))x<Lt  =  [  r(t  -  T)vt[x,r)dT 

Jo  Jo 

=  r(0)u(x,t)  -  r(t)v0(x)  +  f  r'(t  -  r)o(x,  r)dr. 

Jo 


(r) 


(2.2) 


Thus,  for  classical  solutions  (2.1)  is  equivalent  to 


'  wt  =  vx ,  vt  =  <p{w)x  +  T[v\,  x  e  R,  t  >  0, 

<  (2.3) 

w(x, 0)  =  to0(x),  v(z,0)  =  v0(x),  x  €  R, 


where 

J\v\{x,t)  r(0)u(x,t)  -  r(t)v0{x)  +  j  r'(t  -  r)v(x,r)dr.  (2.4) 

Jo 

Under  physically  reasonable  assumptions  (e.g.,  k  =  o',  where  a  :  [0,  oo)  — ►  R  is  smooth, 
positive,  decreasing  and  convex,  r(0)  =  o'(0)  <  0)  the  term  r(0)o(x,f)  has  a  damping 
effect.  Its  effect  is  dominant  close  to  equilibrium  (u;  =  0,v  =  0,0(0)  =  0),  thereby 
inducing  property  (a)  ([22], [11], [39]).  However,  for  suitably  large  data,  the  hyperbolic  part 
of  (2.1)  dominates  irrespective  of  the  sign  of  the  memory,  and  property  (b)  holds  if  0  is 
convex  or  concave  ([8], [31], [35]). 

To  discuss  the  results  for  classical  solutions  of  (VE)  when  0  ^  0,  we  assume  that  the 
smooth  constitutive  functions  0,0  defined  on  R  satisfy 


0(0)  =  0(0)  =  0,  0'(-)  >  0,  0'(O)  >  0. 
4 
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To  simplify  the  exposition,  we  also  assume  that  the  kernel  k  =  a',  where  a  :  ]0,oo)  R  is  a 
smooth,  positive,  nonincreasing,  and  convex  function  on  [0,  oc);  such  kernels  are  relevant  for 
applications  in  viscoelasticity,  (more  generally,  one  can  assume  that  a  is  strongly  positive 
on  f0,oc),  [33], [ll], [12], [16], [39;).  The  following  results  regarding  the  behavior  of  classical 
solutions  of  (VEj  have  been  established  under  additional  technical  assumptions  which  are 
omitted: 

(a)  If  d>'(0)  -  a(0)iy'(0)  >  0  (i.e.,  the  equilibrium  stress  modulus  is  positive)  and  if  the  initial 
data  wo,  Vo  are  sufficiently  smooth  and  small,  it  is  shown  in  [12],  [  16]  using  delicate  energy 
estimates  and  properties  of  Volterra  kernels  that  (VE)  has  a  globally  defined  smooth 
solution  which  decays  as  t  — ►  +oo.  Similar  global  results  for  more  general  models  arc 
discussed  in  [15].  [37;  Ch.  IVj. 

(b)  By  contrast,  for  sufficiently  large,  smooth  data  w0,  v0  and  smooth  kernels  k  irrespective 
of  sign  properties,  the  solutions  of  (VE)  develop  singularities  in  finite  time;  the  first  deriva¬ 
tives  of  w  and  v  become  infinite  while  w,v  remain  bounded.  This  suggests  the  formation 
of  a  shock  front.  Such  results  are  established  for  cj)  convex  or  concave  by  the  method 
of  characteristics  in  [8]  and  [31],  Similar  results  for  other  integral  models  with  smooth 
memory  kernels  and  for  several  differential  models  are  given  in  [37],  Chapter  II.  Numerical 
evidence  of  finite-time  breakdown  of  smooth  solutions  for  large  data  and  for  the  onset  of 
a  shock  structure  for  (VE)  are  provided  in  [26]. 

A  parallel  theory  regarding  properties  (a)  and  (b)  has  been  developed  in  [30],  [25],  [7], 
[23]  for  the  simpler  model  Cauchy  problem  for  the  single  conservation  law  with  memory 

wt  +  ox  —  0,  u;(x,0)  =  w0(x),  ( CLM ) 

where  a  is  given  by  (CA).  When  k  =  0,  (CLM)  reduces  to  the  Burgers  equation 

wt  +  ^(w)*  =  0,  (S) 

while  the  analogue  of  (FE)  is 

wt  -I-  <t>{w)x  +  w  =  0.  ( FB ) 

Our  discussion  of  (VE)  is  limited  to  kernels  which  are  smooth  on  [0,  oo).  However, 
we  remark  that  a  result  confirming  property  (a)  for  (VE)  with  a  singular  kernel  (k  =  a'. 


a  positive,  decreasing,  convex  on  [0,oo),  and  a'  having  an  integrable  singularity  at  zero) 
has  been  established  by  Hrusa  and  Renardy  [17]  (also  see  these  proceedings,  paper  by 
Renardy),  if  i  is  confined  to  a  finite  interval;  it  is  not  known  whether  property  (a)  holds 
for  the  corresponding  Cauchy  problem.  Moreover,  it  is  not  known  whether  property  (b) 
holds  for  (VE)  with  such  singular  kernels. 

The  results  (b)  discussed  above  provide  a  motivation  for  studying  weak  solutions  of 
(VE).  So  far,  only  special  results  concerning  weak  solutions  of  (VE)  have  appeared  in  the 
literature.  Greenberg  [14]  (also  see  [37],  Section  2.6)  established  the  existence  of  travelling 
wave  solutions  (steady  compression  shocks)  for  a  history  value  problem  associated  with 
(VE).  Regarding  unsteady  weak  solutions,  we  present  a  new  existence  result  for  the  Cauchy 
problem  (2.1)  when  the  data  u?o,uo  €  L°°(R),  and  not  necessarily  small. 

We  assume:  the  constitutive  function  <p  satisfies 

{<p  :  R  — >  R  is  a  twice  continuously  differentiable 
function  such  that  <p'(w)  >  0,  w  £  R;  . 

p  has  a  single  inflection  point  at  w  =  W{  and  is  '  ‘  ' 

convex  on  (u/,-,  oo)  and  concave  on  (— oo,u>t). 

The  kernel  k  satisfies 

k  :  [0, oo)  — ►  R,  fee  C'[0,oo),  (2.6) 

and  the  data  tuo^),  v0(:r)  satisfy 

u>0(*),  *o(*)  6  L°°{ R)  n  L2(R).  (2.7) 

The  following  result  establishes  property  (c)  for  the  system  (2.1). 

Theorem  2.1.  Let  the  hypotheses  (2.5)-(2.7)  be  satisfied.  Given  T  >  0,  there  exists  a 
weak  solution  {w(x,t),v[x,t)}  of  (2.1)  on  R  xfO,!1],  such  that 

(w,  v)  €  L°°( (0,  TJ;  L2(R))  n  L°°(R  x  [0,  T]). 

The  proof  is  sketched  in  Section  3  using  the  method  of  compensated  compactness  of 
Murat  [27]  and  Tartar  [40,41.42];  complete  details  of  the  proof  will  appear  in  {32].  This 
approach  has  been  employed  with  success  by  Tartar  [41]  to  obtain  L°°  solutions  for  the 
general,  scalar  Burgers  equation  (B),  by  DiPerna  [13]  and  Rascle  [36]  to  construct  L°° 


a 


«u 


solutions  to  the  hyperbolic  system  (E),  and  by  Dafermos  [10]  to  obtain  L°°  solutions  for 
(CLM);  unfortunately,  the  elegant  approach  of  [10]  does  not  seem  to  apply  to  (VE)  or 
(2.1).  We  also  remark  that  Boldrini  [2]  used  the  method  of  compensated  compactness  to 
show  that  as  the  memory  weakens  (i.e.,  in  (2.1)  k  ~  k(6,t)  =  0(6)),  uniformly  bounded 
solutions  of  (2.1)  converge  to  a  weak  solution  of  (E)  as  6  J.  0. 

Serious  difficulties  remain  to  be  overcome  in  order  to  prove  a  result  similar  to  The¬ 
orem  2.1  for  (VE)  when  rp  ~  <P-  The  central  issue  is  whether  the  memory  term  /J  k(t  - 
t)xP(w)(x,t))  zdr  can  be  shown  to  be  of  lower  order;  such  is  shown  to  be  the  case  if  xp  =  <p, 


even  for  weak  solutions. 


3.  Sketch  of  the  Proof  of  Theorem  2.1. 


The  weak  solutions  of  (2.1),  in  the  space  of  bounded  measurable  functions,  will  be 
constructed  as  limits  of  solutions  of  the  regularized  (parabolic)  system 


Wt  =  Vx  +  £WXX 

Vt  =  <p(w)x  +  fek(t  -  r)<p(w(x,r))xdr 

+e(v*x  +  /0*  k(t  -  t)vxx(x,  r)dr),  (x,  t)  6  R  x  (0,  T] 


w(x,0)  =  w0(x),  v(x,0)  =  v0(x), 


x  g  R, 


as  the  real  parameter  e  j  0.  This  unconventional  regularization  preserves  one  of  the  main 
features  of  (2.1);  namely,  that  the  integration  with  respect  to  t  offsets  differentiation  with 
respect  to  r  and  the  memory  terms  are  of  lower  order.  Its  relevance  is  revealed  by  the 
following  calculation.  Convolving  (3.1) 2  with  r(<),  and  using  equations  (r)  and  (2.4),  we 


/  k(t  -  T)(ip(w(x,r))x  +  evzx(x,T))dT  =  f  r(t  -  r)vt(x,T)dT 

Jo  Jo 


=  r(0)v(x,t)  -  r(t)v0(x)  +  [  r'(t  -  r)v(x,  r)dr  =  7\v\(x,t). 

Jo 


Thus,  the  initial  value  problem  (3.1)  can  be  written  in  the  form 


=  v*  +  ewzx,  =  <p[w)x  +  7\v)  +  evxx,  ( x,t )  G  R  x  (0,T] 


;(x,0)  =  tuo(i),  t'(x,0)  =  v0(x),  x  G  R. 


It  is  evident  that  systems  (3.1)  and  (3.3)  are  equivalent  for  smooth  data  and  classical 
solutions.  For  the  data  in  LC>D(R)  n  L2(R),  the  above  ideas  extend  to  weak  solutions; 
i.e.,  it  is  showm  in  32]  that  (3.1)  and  (3.3)  (also  (2.1)  and  (2.3))  are  equivalent  for  L°° 
solutions,  using  Cutest  functions  which  have  compact  support  on  R  x  [0,  T]  for  any  T  >  0 
and  which  vanish  at  t  =  T.  Moreover,  (2.2)  and  (3.2)  hold  in  the  weak  sense. 

Our  next  objective  is  to  obtain  a-priori  estimates,  independent  of  e,  for  solutions 
{w£  ,v£}  of  (3.1)  corresponding  to  initial  data  w0,  t>o  6  L°°(R)  n  L2(R);  we  will  employ  the 
equivalent  system  (3.3).  The  following  result  is  needed  to  complete  the  proof  of  Theorem 
2.1. 

Theorem  3.1.  Under  the  hypotheses  (2.5)-(2. 7),  for  each  e  >  0,7'  >  0,  the  initial 
value  problem  (S.S)  (respectively  (S.l))  has  a  unique  solution  {w£(x,t),v£(x,t)}  defined  on 
R  x  [0,T]  such  that  w£,v£  G  C([0,  T];  L2(R))  n  L°°(R  x  [0,T]),  w\,  v£,  wxx,  vexx,  w\,  v£  G 
C((0,  T}\ L2(R))  and,  also,  u)x,vx  G  L2(R  x  [0,  T]).  In  addition,  as  e  [  0,  the  families 
{w£(x,t),ve(x,t)}i:>o  and  {ex^2w£x{x,t),£l^2vx{x,t)}e>o  lie  in  bounded  sets  of  L°°([0,T]; 
L2(R)  n  L°°(R))  and  L2( R  x  [0,  Tj),  respectively. 

The  existence,  uniqueness  and  regularity  properties  of  the  solution  ( w£ ,  v£)  of  (3.3)  are 
established  in  [32]  by  a  fairly  standard  fixed-point  argument  on  an  appropriately  chosen 
Banach  space  using  energy  methods  and  the  L2-theory  of  the  heat  equation.  We  will 
sketch  the  proof  of  the  a-priori  estimates;  we  emphasize  that  due  to  the  nonlocal  nature  of 
the  memory-term  in  (3.3),  the  L°°-estimates  cannot  be  obtained  by  finding  an  invariant 
region  [4]  (compare  with  the  proof  of  the  existence  of  weak  solutions  for  (E)  in  [13]).  To 
avoid  a  confusing  notation,  we  now  drop  the  superscript  e. 

The  a-priori  estimates  are  deduced  by  using  a  class  of  exponentially  growing  convex 
entropies  for  (E)  constructed  by  Dafermos  [9].  Let  {i^(i,  t),  v(x,  f)}  be  the  solution  of 
(3.3)  on  R  x  [0,  T]  satisfying  the  regularity  properties  of  Theorem  3.1;  this  assumption 
justifies  the  calculations  which  follow.  In  the  sequel,  C  will  stand  for  a  generic  constant 
depending  on  the  L°°(R)  and  L2(R)-norms  of  the  initial  data,  on  the  C^O, Tj-norm  of 
r(t),  on  properties  of  the  function  <p{w),  on  T,  but  not  on  e.  Whenever  the  constant 
depends  on  e  it  will  be  denoted  by  Cc.  The  analysis  of  (3.3)  will  be  based  on  the  concept 
of  entropy-entropy  flux  pairs  for  the  elastic  problem  (E)  (cf.  Lax  [20]).  A  smooth,  convex 


function  rjlw,  u)  defined  on  RxR  is  an  entropy  for  (E),  with  corresponding  entropy  flux 
q{ w,v),  if 

dtr](v)(i,t),v(i,t))  +  dzq(w(x,t),v(x,t))  -0  (3.4) 

for  any  smooth  solution  {te(x.  t),  r(x,  t)}  of  (E).  Such  pairs  are  generated  as  solutions  of 
the  system  of  equations 

in  - 

(3.5) 


|  Qw  -  ~<p'{w)riv 
\  Qv  =  — ^?UM 

provided  77(11/,  u)  is  convex.  Eliminating  q[w,v)  in  (3.5),  we  find  that  rj(u>,v)  must  be  a 
convex  solution  of  the  linear  wave  equation 


uiw  —  P  {,w}t)vv\  (3.6) 

q(w,v)  is  then  determined  by  (3.5).  A  classical  example  of  an  entropy-entropy  flux  pair  is 

1  fw 

*?(«>, v)  =  -v2  4-  J  p(£)d$,  q(w,v)  =  -vp(w).  (3.7) 

For  each  entropy-entropy  flux  pair  {77(11;,  v),  q(u>,  t/)}  for  (E),  we  denote  by  (17(11;,  t), 
9(ui,t’)}  the  parts  of  the  pair  {r),q}  which  vanish  to  quadratic  order  in  (u;,t;).  Then  {17 ,9} 
also  form  an  entropy-entropy  flux  pair  for  (E),  and  a  simple  computation  yields  the  identity 
dtfj{w,v)  +  dxq(w,v)  =  fjv(w,v)?\v]  +  ed2zf](w,v) 

-  v)w 2  +  2r]wv(w,  v)wxvx  +  tjvv{w,  v)vl}. 

Integrating  (3.8)  over  Rx(0,  t),  0  <  t  <  T,  we  obtain 


(3.8) 


/: 


T](w(x,t),v(x,t))dx 

5 

•  t  r  OO 


+  £  /  /  (riww(w,v)wl  +  2rjwv(w,v)wxvx  +  r}vv(w,v)vl)dxdr 

Jo  J  —  oo 

/oo  r  t  f  OO 

77(u/0(x),  u0(x))dx  +  /  /  T7„(u,*,  v)T\v)dxdr,  0  <t<T. 

-00  J 0  J  —  00 


(3.9) 


For  the  first  estimate,  we  employ  the  pair  (3.7)  and  use  (3.9)  to  deduce 

/OO  ,  ft  fOO 

[~v2(x,t)  -r  $(u/(x,  t))]dx  -i-  e  /  /  (<p'[w)w2  +  v2)dxdr 

■  00  “  J  0  J  -  00 

=  f  (J^(x,0  +  $(u/o(x,0)jrf* 

7—00  ^ 
rt  roo 


+  f  j  v?\v]dxdT ,  0  <  t  <  T, 
Jo  J — 00 


(3.10) 


v*  ere 


$(w)  :=  [  (<p(£)  -  vp(Q))^. 
Jo 


Assumption  (2.5)  implies  that 


<p'(u>,)  2 


<p'(ui)  >  <f'(wi)  >  0,  and  >  — - — -  w 2  >  0,  for  every  w  €  R. 

mt 

Using  (2.4)  and  the  Cauchy-Schwarz  inequality,  a  simple  calculation  yields 


Combining  (3.10)  with  (3.12)  and  (3.13),  gives 


/oo  ft  roo 

u2(x, t)dx  <  C  +  C  /  I  v2dxdr, 
-OO  J  0  J  —  oo 


whence  by  Gronwall’s  inequality,  (3.10)  and  (3.12), 


/oo 

[u2(x,t)  +  w2(x,t))dx 

-oo 

4-  £  f  I  [w2  +  v2x)dxdr  <  C,  0  <  t  <  T, 
Jo  J- oo 


(3.11) 


(3.12) 


/  f  v7[v]dxdr  <  C(  1  +  f  f  v2(x,  r)dxdr),  0  <  t  <  T.  (3.13) 

Jo  J ~ oo  Jo  J  —  oo. 


(3.14) 


(3.15) 


where  C  is  independent  of  e.  Thus  {u;(x,t),t;(x,t)}  lies  in  a  bounded  set  of  L°°{\ 0,T]; 
L2(R))  and  {el^2wx(x,t),e1^2vx(x,t)}  lies  in  a  bounded  set  of  L2(R  x  [0,r]),  independent 


a 

•Sa 

m 


The  proof  of  the  Z/°°-estimates  independent  of  e  is  more  subtle.  For  this  purpose,  v 
extend  the  development  of  Dafermos  [9j.  The  following  facts  are  proved  in  [9,  Section  2]: 
For  <p{w)  as  in  (2.5),  the  wave  equation  (3.6)  admits  a  class  of  solutions  {^(^(u;,  u)}fc>o  on 
RxR  which  are  strictly  convex  and  grow  exponentially  at  infinity.  These  solutions  have 


the  form 


Tj(k\w,v)  =  Y^[w)cosh  kv ,  0  <  k  <  oo, 


where  Y^(w)  is  the  solution  of  the  initial  value  problem 

f  Y^k^'{w)  =  k2ip'(w)Y^(w), 
U^^l,  F(fc)'N=0,  0  <  k  <  oo. 


(3.16) 


(3.17) 


iJS 

SR 

I 

:i'N 


The  functions  y^’c,(t u)  satisfy  the  estimates 

y(fc)(u>)  >  cos/i;A:y/^,(u,’,)(uj  -  u.’,)],  -oo  <  u;  <  oc,  0  <  k  <  oo,  (3.18) 

\Y[k)  (id) |  <  k\/<p'(w)Y ^(tc),  -oc  <  w  <  oo,  0  <  k  <  oo  (3.19) 

and 

ru  _ 

Y^k\w)  <  exp[k  /  \/'-p'{£)d£]i  -oo  <  tu  <  oo,  0  <  <  oo.  (3.20) 

J  Ult 

We  will  estimate  the  solution  (u.(x,t),u(x,f)}  of  (3.3)  by  monitoring  the  evolution  of 

/oo 

^C/c)(lD(x,f),u(x,t))dx. 

-oo 

In  view  of  the  convexity  of  v),  (3.9)  yields 

/oo  r  oo 

fj(*)(w(x,t),v(x,t))dx  <  I  r?^(u;o(x),uo(x))dx 

ft  ,00  .  (3-21) 

+  /  /  i?|,fc)(tD(x,r),t;(x,r))/{v](x,r)dxdr,  0  <t<T, 

Jo  J- 00 

where,  by  (3.16), 

fjW{w,v)  =  V^{w,v)  -  y  ^(0)  -  y(fc)'(0)u;  >  0,  (3.22) 


rjlk\w,v)  —  k{tanh  kv)rj^k^(w,v). 


(3.22) 


(3.23) 


Using  (3.23),  (3.22),  (3.19),  the  Cauchy-Schwarz  inequality  and  (2.4),  a  lengthy  calculation 
shows  that  the  last  integral  in  (3.21)  can  be  estimated  as  follows: 

|/  [  f)[k](xv,v)J\v)dxd,T 

1 Jo  J-00 

noo 

k(tanh  kv)(fj^(w,v)  ■+■  y^fc^(0) 

-00 

+  y(fc)  (0)u>)  J\v]dxd,T  J 

<  CkWk >(0)  J‘  J°j\w\  +  !»|)|/[»]| dxir  (3.24) 

+  k  f  f  fj^(u;, v)| 7\v]\dxdr 

Jo  J-OC 

ftfO o  f  00 

<Ck2Y^k\0)[l  I  (u>2  +  v2)dxdr  +  /  u2(x)dx] 

J 0  J -00  J —00 

+  Ck  f  M{r)[f  fj(kS>(w(x,T),v(x,T))dx]dT, 

Jo  J-00 


m 


where 


luri 


c 


'3 

it 


v>  UCIC 

M(0  :=  sup {1  +  jt/(x,f)J  +  /  1  v[x,T)]dr},  0  <  t  <  T.  (3.25) 

rGR  Jo 

Combining  (3.21),  (3.24),  (3.15),  and  then  using  the  Gronwall  equality  yields  the  estimate 


f  rj(k)(w(x,t),v(x,t))dx 

J  —  OO 

<{[  fjW(w0(x),vo(x))dx 

+  Ck2Yw[0)}exp{Ck  [  M(f)dr},  0  <  t  <  T. 

Jo 


(3.26) 


Raising  both  sides  of  (3.26)  to  the  l/k  power,  letting  k  — *  oo  and  then  taking  the  logarithm 
of  both  sides  of  the  resulting  inequality  (see  [32]  for  justification  of  this  procedure),  we 


conclude 


sup{l  -f  |tn(x,t)|  +  |v(x,t)|}  <  C  +  C  [  M(r)dr,  0  <t<T.  (3.27) 

x£R  Jo 


Substituting  (3.25)  into  (3.27)  yields 


S(f)  <  C  +  C(l  +  T)  /*  S(r)dT, 

Jo 


(3.28) 


where  S(t)  :=  sup{l  +  |u;(x,t)|  +  |v(z,t)]}.  Finally,  integrating  (3.28)  we  conclude 


sup{l  -f  \w(x,  t)l  +  \v(x,t)\}  <  Af,  0  <  t  <  T, 
z€R 


(3.29) 


where  the  constant  Af  depends  on  T  but  not  on  e.  Therefore,  the  solution  {u>(x,t),r(x,t)} 
of  (3.3)  is  bounded  on  R  x  [0,7’],  uniformly  in  £  >  0.  This  completes  the  sketch  of  the 
proof  of  the  a-priori  estimates  in  Theorem  3.1.  ■ 

Our  final  task  in  proving  Theorem  2.1  is  to  use  Theorem  3.1  to  justify  passage  to 
the  limit  as  e  j  0+  in  the  Cauchy  problem  (3.1).  For  e  >  0,  let  {u>e(x,t), ve(x,t)}  be 
the  solution  of  (3.1)  on  Qt  :=  R  x  [0,7’],  with  regularity  properties  as  in  Theorem  3.1. 
By  virtue  of  (3.29)  the  family  of  functions  {iue(x,  t),  ve(x,  £)}e>o  is  uniformly  bounded  in 
£°°(Qr)-  There  exist  functions  w(x,  t),  v(x,  t)  and  <p(x,t)  in  L°°(Qt )  such  that,  along  a 
subsequence,  w£(x,t)  w(x,  t),  ve(x,t)  v(x,t),  and  v?(t ue(x,t))  <p(x,t)  in  i°°-weak 
star  as  e  J.  0.  The  objective  is  to  show  that  {u;(x,  t),  v(x,  <)}  is  a  solution  of  (2.1)  in  the 
sense  of  distributions  for  (x,t)  €  Qt- 


The  major  obstacle  to  overcome  is  that,  in  general,  nonlinear  functions  are  not  con¬ 
tinuous  under  weak  star  convergence,  and  it  does  not  follow  that  <p(x,t)  =  tp(xv(x,t)). 
Under  weak-star  convergence,  such  composite  weak  limits  have  been  characterized  by  Tar¬ 
tar  [40,41,42.  as  expected  values  of  a  family  of  probability  measures  V(x,t)i  called  the  Young 
measures.  The  relevant  issue  here  is  whether  the  Young  measure  reduces  to  a  Dirac  mass. 
Using  the  theory  of  compensated  compactness  and  a  class  of  entropy-entropy  flux  pairs  in¬ 
troduced  by  Lax  [20:,  DiPerna  [13;  proved  the  following  key  result  which  will  be  employed 
in  the  sequel. 

Proposition  3.2  (DiPerna).  Let  {iy£,u£}  :  Qt  —>  R  be  a  collection  of  functions  such  that 

llu;4r|U00(Qr)  ^  <  C, 

where  C  is  a  constant  independent  of  e.  Suppose  also  that  for  any  smooth  entropy- entropy 
flux  pair  {q(u>,  r),  q(w,  r)}  of  (E)  with  <p  satisfying  (2.5), 

dtrf(we(x,t),vc(x,t))  +  dxq(we(x,t),ve(x,t)) 

lies  in  a  compact  set  of  HfcliQr)-  Then  there  exists  a  subsequence  {wc  ,ve  }  and  functions 
w,v  £  L°°(Qt )  such  that 

wc  [x,t)  — ►  w(x,t),  v€  (x,t)  — ►  v(x,t),  a.e.  for  (x,t)  £  QT 

as  e'  i  0. 

We  apply  Proposition  3.2  to  the  family  {w€(x,  t),  v£(x,  t)}£>o  of  solutions  of  (3.j) 
(which  also  satisfy  (3.3)).  A  straightforward  calculation  using  (3.3)  and  (3.5)  yields 

dtr)(we,  vc)  +  dxq(we,v£) 

=  £1/2dx(e1/2riw(wc ,  ve)wex  +  el'2r,v(we ,*>') 

(3.30) 

-  e[riww{we  ,ve){wx)2  +  2r)wv(wc,ve)wezvex  +  rjvv(w£,  ve)(vx)2] 

+  r)v(we,ve) 7{vc\  :=  /i  -  /2  +  J3> 

The  a-priori  estimates  (3.15)  and  (3.29)  imply  that  the  family  I\  converges  to  0  and  is 
thereby  compact  in  H~1[Qt),  while,  the  family  I 2  is  bounded  in  L1{Qt)-  In  addition, 


lies  in  a  bounded  set  of  IV  1,00(Qr)-  Using  a  lemma  of  Murat  [28]  we  deduce  that  (3.31) 
lies  in  a  compact  set  of  Hlol(Qr )•  Finally,  since  /3  1  ies  in  a  bounded  set  of  L2(Qt)  and 
thus  in  a  compact  set  of  H[~0'(Qt),  the  left  side  of  (3.30)  also  lies  in  a  compact  set  of 
//,;<! (Qr)-  Thus,  Proposition  3.2  implies  that 

we  (x,t)  — ►  u>(x,t),  v£  ( x,t )  — *  v(x,t)  a.e.  in  Qt, 

along  a  subsequence  s'  j  0,  and  permits  passage  to  the  limit  e  j.  0  in  (3.1)  in  the  sense 
of  distributions.  The  pair  of  functions  {u>(x,J),r(z,*)}  belongs  to  L°°('0,  T];  L2(R))  n 
L°°(Qt )  and  is  a  weak  solution  of  (2.1)  in  Qt-  This  completes  the  sketch  of  the  proof  of 
Theorem  2.1.  ■ 
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